Dissipative quantum phase transition in a quantum dot 
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We study the transport properties of a quantum dot (QD) with highly resistive gate electrodes, 
and show that the QD displays a quantum phase transition analogous to the famous dissipative 
phase transition first identified by S. Chakravarty [Phys. Rev. Lett. 49, 681-684 (1982)]; for a 
review see [A. J. Leggett et at, Rev. Mod. Phys. 59, 1 (1987)]. At temperature T = 0, the charge 
on the central island of a conventional QD changes smoothly as a function of gate voltage, due 
to quantum fiuctuations. However, for sufficiently large gate resistance charge fluctuations on the 
island can freeze out even at the degeneracy point, causing the charge on the island to change in 
sharp steps as a function of gate voltage. For Rg < Rc the steps remain smeared out by quantum 
fiuctuations. The Coulomb blockade peaks in conductance display anomalous scaling at intermediate 
temperatures, and at very low temperatures a sharp step develops in the QD conductance. 



The single electron transistor (SET) is one of the most 
basic mesoscopic devices: A conducting island or quan- 
tum dot is attached by tunnel barriers to two leads and 
a capacitively-coupled gate electrode sets the number 
of electrons on the dot. For low enough temperatures, 
T <C Ec, charge fluctuations of the dot are suppressed 
except when the gate is tuned to make two charge states 
nearly degenerate. At these "charge degeneracy points" 
the charge on the dot strongly fluctuates. For typical 
metallic SETs with a very large number of tunneling 
modes quantum fluctuations of the charge turn out to 
be suppressed at low temperatures Xi ^ M- For semi- 
conducting SETs with single mode junctions, however, 
quantum fluctuations of the charge are important and 
broaden out the charging steps at low T: In the limit of 
vanishing level spacing, — > charge fluctuations are 
described by the two-channel Kondo model , while for 
T <C (Se one recovers the so-called "mixed valence" regime 
of the Anderson model [^. 

In the above discussion we neglected the effect of 
Ohmic dissipation in the lead electrodes. While this 
has been extensively studied for SETs with a very large 
number of tunneling modes Q , there is much less known 
about the effects of dissipation in the Kondo regime: In 
a recent paper Le Hur showed that, assuming a contin- 
uum of quantum levels on the SET and a single tunnel 
mode, coupling to a dissipative bath drastically modifles 
the results of Ref . 0| : large enough dissipation drives a 
Kosterlitz-Thouless-type phase transition and leads to a 
complete suppression of charge fluctuations even at the 
degeneracy point 0,Q|- However, for most semiconduct- 
ing devices the level spacing 5e cannot be neglected in 
comparison to temperature, and spin fluctuations must 
also be considered. Here we shall therefore investigate 
the effects of dissipation at temperatures far below the 
level spacing on the dot, T ^ 5e, a more realistic low- 
temperature limit for typical semiconductor SETs. As we 



show below, a dissipation-induced quantum phase tran- 
sition takes place for T ^ (5e as well, although with dif- 
ferent and more complicated properties due to the in- 
terplay of charge and spin fluctuations, and at a larger 
dissipation strength (gate resistance) than that needed 
for <5e->0 Aia]- 

In this T <C (5e regime, the coupling of the quantum 
dot to the gate voltage is usually described by the Hamil- 
tonian. 
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where Eq = e^/2Ci; denotes the charging energy, with 
Cs the total capacitance of the dot, and e the electron 
charge. We retain only one single-particle level d, and 
we assume that it is empty or singly occupied, depend- 
ing on the dimensionless gate voltage, Ug y. Assuming 
weak coupling between the dot and the source and drain 
electrodes, charge transfer can be described within the 
tunneling approximation. 



(2) 



were ipa- (e) annihilates an appropriate linear combination 
of left and right lead electrons of energy e that hybridize 
with the dot state do-, and satisflcs the anticommutation 
relation {Me),iP^,{e')} = S^^^Sie-e') 0- 

Throughout 

this paper we assume that the quantum dot is close to 
symmetrical but our analysis carries over easily to asym- 
metrical dots as well. 

Eqs.O and 10) are thought to provide a satisfactory 
description of the SET for T <ti 6e for most experimen- 
tal situations studied so far, including the Kondo regime 
0. However, Eq. |^ does not account for the relax- 
ation of electrostatic charges in the nearby electrodes: in 
reality, when an electron tunnels into the dot, an elec- 
trostatic charge SQ = eCg/Cj: is also generated on the 
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gate. Transferring this charge from the outside world to 
the gate electrode through a shunt resistor requires time, 
and creates dissipation llj. Consequently, tunneling be- 
tween dot and leads will be suppressed by Anderson's 
orthogonality catastrophe. The simplest way to account 
for this shunt resistance is to add a term Q 



diss 



= A 



(3) 



where the operator n = dl-d^ measures the number 
of electrons on the dot, and the bosonic field ip describes 
charge density excitations in the gate electrode partici- 
pating in the screening process, with their imaginary time 
correlation function defined as (Tr(^(r)(^(0)) = [l^ . 
The dissipation strength a = AV 2 can be estimated fol- 
lowing the procedure of Ref. |l3j , and we find 
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with Rg the low frequency resistance of the gate elec- 
trode, and Rq = h/2e'^ the resistance quantum. In the 
absence of dissipation, Eqs.QI and ^ would predict that 
the average charge on the island changes smoothly from 
(n) « to {fi) « 1 around Ug w 1/2 due to the presence 
of strong quantum fluctuations at ng w 1/2 5]. As we 
shall see below, the properties of this transition may be 
dramatically modified for a > 1/2, producing a sharp 
step in (h) at T = 0. 

Since we shall focus on the vicinity of the charging step, 
rig « 1/2, we restrict the Hilbert space to the empty state 
|0) and the two singly-occupied states \a) = (ij.|0) of the 
dot, and approximate the sum of Eqs. |5J) and (jSJ as 



(5) 



where A = Ec{l — 2ng) measures the difference between 
the classical energy of the two charge states of the dot, 
the operator Q= (| T)(T l + l iXi | - 2|0)(0|)/3 measures 
the difference between the charge on the dot and its value 
at the charge degeneracy point, and ipa — J V-'crC^)- 

To map out the complete phase diagram of the model, 
we bosonize the above Hamiltonian 14], and perform a 
renormalization group analysis using an operator prod- 
uct expansion method Although we use different 
methods, the derivation of the scaling equations is simi- 
lar to that in Ref. 16]. In fact, a mapping between the 
present model and the one studied in Ref. enables 
us to carry over much of the analysis as well. Intrigu- 
ingly, a new exchange coupling j is also generated under 
scaling [l^: 
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where the spin operators denote = \ Yla p \^)^'ai3{P\^ 
with cr' the Pauli matrices. To lowest order, the scaling 



equations read 0] 
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where Z = In a is the scaling variable, and v — Va^^"^ 
is the dimensionless tunneling rate with a a microscopic 
time scale initially of the order of ao ~ l/5t. The cou- 
pling A = aA is the dimensionless splitting of the charge 
states, slightly renormalized by the dissipative coupling. 
Note that these equations are exact in the dissipative 
coupling a, but they contain contributions in the small 
couplings V and j only up to second order. 

Following Ref. 0] , we first neglect the effect of the ex- 
change coupling j and focus on the vicinity of the degen- 
eracy point, defined through the condition A(Z — > oo) — 
0. As we shall see, this approximation describes much of 
the experimentally-relevant regions, because the Kondo 
temperature associated with spin fluctuations turns out 
to be usually quite small in the vicinity of the phase 
transition, a ~ ac, even for single occupancy. Within 
this approximation we can solve the scaling equations 
analytically 18], and for large enough a we find that as 
T — > the effective tunneling v scales to zero, while a 
scales to a finite value aoo larger than 1/2. Thus, quan- 
tum fluctuations of the dot charge vanish as T 0. In 
the above approximation, the relation between the cor- 
responding critical value of the parameter a = ai"^ and 
the level width F « 2TTv'^Se is determined by the equation 
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and is plotted in Fig. ^ For small values of v we find 
a''c^\v) = 1/2-1- + v'^ + . . . , and the transition is 
of Kosterlitz-Thouless type the localized side. 



a > a. 



(0) 



(or F < Fc""*), at the degeneracy point the 
height 5G of Coulomb blockade peaks scales to zero as a 
power law |1S] . 
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with Gq the quantum conductance. On the metallic 
side, on the other hand, quantum fluctuations always 
dominate and preserve conductance even at T — 0, 
though near the transition the conductance shows a non- 
monotonic behavior: 5G{T) first slowly decays and then 
starts to increase below a temperature T* that vanishes 
exponentially as one approaches the phase transition. 
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FIG. 1: Schematic phase diagram of the SET in the presence 
of dissipative coupling, denotes the critical value of a, 
while a'"' is its value obtained by neglecting the generated 
exchange coupling j. For a > Oc there is a phase transition 
from (n) = to (n) = 1, while in the more familiar situation 
of weak dissipation there is a crossover. 



T* « 5e exp{-7r/2(a^°^ - a^)^/^}, until finally a mixed 
valence state with a large conductance is formed at a 
temperature T** ~ T*"^ /5t < T* . For the critical value 
of a, SG decays to zero logarithmically, 
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(0)^ 



G, 



1 



Q 



\n\6e/T) 



(13) 



The gate- voltage dependence of the expectation value 
of the charge on the dot, (Q), can be directly detected by 
another SET or point contact electrostatically coupled 
to the dot. So far, we only considered the degeneracy 
point, A(oo) = 0. However, in the localized phase, A is 
a relevant perturbation and generates a first order phase 
transition at T = (see Fig. QJ. Correspondingly, (Q) 
jumps [Tsf at the degeneracy point. At finite tempera- 
tures the width Srig of this charge step becomes finite, 
5ng ~ T/Ec. 

On the delocalized side v scales to large values and 
a — > so that the dissipation effectively disappears from 
the problem, and we recover the well-understood mixed 
valence fixed point of the Anderson model [^H^. At this 
fixed point the charge susceptibility is finite, and thus the 
charging step remains smeared out by quantum fluctua- 
tions even at T = 0. The width of the step is roughly de- 
termined by the scale T**, which vanishes exponentially 
fast as we approach the transition. The corresponding 
T ~ phase diagram is sketched in Fig. 

We shall now discuss how the thus-far neglected ex- 
change coupling j changes the above picture. Let us 
again first examine the degeneracy point, A — + 0, and let 
us assume that we are deep in the localized phase. Here 
charge fiuctuations (v) are irrelevant, however a small ex- 
change coupling is still generated at the beginning of scal- 
ing, and produces a spin Kondo effect at the degeneracy 
point at some Kondo temperature Tk- The scale Tr- is 



typically small, but could likely be pushed into the mea- 
surable range by gradually opening up the quantum dot. 
Although the exchange coupling somewhat renormalizes 
f , it will typically not make it relevant. In other words, 
in the localized phase, at the degeneracy point one scales 
at low T to a special quantum state with lar ge spin fluc- 
tuations but suppressed charge fluctuations [la>l21|- As 
a consequence, at the degeneracy point the conductance 
of the SET becomes of the order of Gq at temperatures 
T ^ Tff, even though charge fluctuations of the SET 
are suppressed and v scales to zero. Despite the strong 
spin fluctuations, A remains a relevant perturbation, and 
therefore, at T = there is still a jump in {n){ng). In- 
terestingly, on the (h) w 1 side of this jump, the Kondo 
effect always develops, the spin is screened, and a Fermi 
liquid is formed for T < Therefore the T = 

conductance is close to the quantum conductance. On 
the other side of the transition, however, (n) « 0, spin 
fiuctuations are suppressed, and no Kondo effect takes 
place. Since the conductance is directly related to the 
scattering phase shifts and thus the occupation number 
(n) through the Friedel sum rule, the T = conductance 
in this localized phase must also have a jump as a function 
of Ug. This regime may be difficult to explore experimen- 
tally since we expect Tk to be small unless both a and 
the tunnel coupling V are sufficiently large. 

Exchange fluctuations play another important role too: 
they increase the value of v and, as a consequence, the 
true critical value of a is somewhat renormalized, a^'^ — > 
Q!c. It is very hard to estimate this critical value of ac, 
but for vanishingly small values of v we have been able 
to show that adv ^ 0) = ai^\v ^ 0) = 1/2. For non- 
vanishing values of v, counter-intuitively, ac seems to be 
shifted to somewhat smaller values |2^ . 

To obtain a more quantitative picture of the phase 
transition, we performed numerical renormalization 
group calculations. We used the Anderson Hamiltonian, 
Eqs. (0 and lO, and represented the field (p by fermionic 
density fluctuations 0- The computed T = charging 
steps are shown in Fig. |21 The steps get sharper and 
sharper as we approach ac, and finally a sudden step 
appears for a > ac- 

The T — temperature AC conductivity is shown in 
Fig.Elin the charge localized phase. The AC conductivity 
clearly shows the Kondo resonance, whose width remains 
finite as one approaches the charge step from the (n) w 1 
side, while a dip of width ^ |A| develops on the (n) w 
side, consistent with the G{uj = 0) conductance having a 
jump. Similar non-monotonic behavior should appear as 
a function of temperature. 

To achieve a large dissipation strength we propose to 
make a shallow two-dimensional electron gas (2DEG) and 
fabricate a highly resistive metallic top gate electrode 
just above the dot, with a resistance larger than Rq. 
According to our estimates, dissipation strength of the 
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peratures an anomalous scaling of the Coulomb blockade 
peaks is predicted. We estimate that this quantum phase 
transition can be detected by coupling a highly resistive 
gate electrode to a SET in a shallow 2DEG. 

We are grateful to P. Simon, K. Le Hur, Q. Si, O. 
Sauret, A. Zaikin and Y. Nazarov for valuable discus- 
sions. This research has been supported by NSF-MTA- 
OTKA Grant No. INT-0130446, Hungarian Grants Nos. 
T046303, NF061726, D048665, and T048782, the Euro- 
pean 'Spintronics' RTN HPRN-CT-2002-00302, and at 
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FIG. 2; Charging steps, computed using NRG for a relatively 
small hybridization. For these parameters Tx/Se ~ 10^^" 
at a « ac. Inset: Temperature dependence of the occu- 
pation number (n) at the critical dissipation, a ~ a^, for 
A/5e = -0.0028800, -0.0028300, -0.0028200, -0.0028198, 
-0.0028190, -0.0028140, -0.0028000 (bottom to top). 



order of a « 1 can be reached in this way. The SET 
can be then tuned through the quantum phase transition 
by either continuously changing the tunneling V ^ or by 
depleting a second 2DEG positioned below the dot and 
thereby changing the total capacitance of the dot and 
hence the value of a. 

In summary, we have shown that sufficiently strong dis- 
sipation in the gate electrodes can drive the SET through 
a quantum phase transition into a state where charge de- 
grees of freedom become localized while spin fluctuations 
lead to a Kondo effect. In this state both the conductance 
and the expectation value of the charge on the SET dis- 
play a jump at temperature T = 0, while at higher tem- 
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FIG. 3; r = AC conductivity of the SET in the localized 
phase for a = 0.75 and T/5t = 0.5 (Go ~ Gq). 
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